LPT-11-102 



Gauge-invariant perturbations at 
second order in two-field inflation 



Eleftheria Tzavara and Bartjan van Tent 

Laboratoire de Physique Theorique, Universite Paris-Sud 11 and CNRS, 
Batiment 210, 91405 Orsay Cedex, France 

E-mail: Eleftheria.Tzavara@th.u-psud.fr, Bartjan.Van-Tent@th.u-psud.fr 

Abstract. We study the second-order gauge-invariant adiabatic and isocurvature pertur- 
bations in terms of the scalar fields present during inflation, along with the related fully 
non-linear space gradient of these quantities. We discuss the relation with other perturba- 
tion quantities defined in the literature. We also construct the exact cubic action of the 
second-order perturbations (beyond any slow-roll or super-horizon approximations and in- 
cluding tensor perturbations), both in the uniform energy-density gauge and the flat gauge 
in order to settle various gauge-related issues. We thus provide the tool to calculate the exact 
non-Gaussianity beyond slow-roll and at any scale. 
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1 Introduction 

The concept of inflationary curvature perturbations was first invoked in order to explain 
the primordial fluctuations that source the CMB anisotropy and structure formation [1-4]. 
The inflationary paradigm has been observationally tested for more than 10 years and its 
prediction for an almost scale invariant spectrum of the first-order curvature perturbations 
has been verified by the data of many experiments (see for example [5]). Since the definition 
of perturbations depends in general on the gauge choice, a gauge-invariant definition of the 
cosmological perturbations is of vital importance to make contact with physical observables, 
which are obviously gauge-invariant. That was investigated in detail in [6] and later in [7]. 
In the mean time the need for more observational quantities than just those based on linear 
perturbation theory has become clear, in order to break the degeneracy of the immense 
number of inflationary models. One of the most fruitful has proven to be the non-Gaussian 
characteristics of the perturbations. This has led to the development of new methods to study 
the combination of the scalar field and metric perturbations, the one sourcing the other, in 
a gauge-invariant way beyond first order. 

It was not until 2003 that Malik and Wands in [8] defined the gauge-invariant quantity at 
second order that reduces to the curvature perturbation in the uniform energy-density gauge. 
In [9] the super-horizon equations of motion of these quantities were derived (but see also 
[10] for a gauge-ready formulation of the perturbations and their equations). Another way 
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to deal with perturbations at second order is the space gradients approach first introduced 
in [11] and specifically the gradient of the fully non-linear curvature perturbation used in 
[12] and defined by the same authors in [13]. The advantage of the method is that when the 
space gradients are expanded to first order they are automatically gauge-invariant. Some 
years later the gradient of the curvature perturbation was redefined in a covariant way in 
[14]. In this paper it was shown that when expanded to second order, this quantity reduces to 
the gauge-invariant curvature perturbation defined in [8] plus a gauge transformation term. 

In this paper we generalize the definition of the gauge-invariant curvature perturbation 
(or the gradient of the relevant fully non-linear quantity) in terms of the energy density to 
a definition in terms of the scalar fields present during inflation and study the consequences 
of this change at second order. Since the scalar fields are the principal quantities during 
inflation, it makes more sense to use these as a starting point, especially in the case of 
multiple-field inflation. Our original need for such a definition was to find the horizon crossing 
contribution to the second-order curvature perturbation in terms of the first-order ones in 
the long- wavelength formalism [15]. Indeed, such a definition helps to fully understand and 
potentially generalize the two formalisms used to compute non-Gaussianity during inflation, 
i.e. the long- wavelength formalism [12, 15, 16] and the 6N formalism [17-21], where instead 
of the energy density, the values of the fields themselves are used. 

In the case of multiple-field inflation not only an adiabatic curvature perturbation is 
produced, but also one or more isocurvature perturbations. We would like to provide the 
same type of study for the isocurvature perturbation as for the curvature one, using the 
generalized quantity deflned in [12], and deduce from that the second-order gauge-invariant 
analogue. This definition of the isocurvature perturbation makes direct contact with the 
scalar fields during inflation (instead of using their pressure) , which we find more useful during 
the period of scalar field domination of the universe. It has a simple physical meaning, that 
is the combination of the fields that remains orthogonal to the field trajectory, as opposed 
to the adiabatic perturbation that is parallel to the field trajectory (and proportional to the 
energy density). 

On a related subject, Maldacena in [22] found the third-order action for the first-order 
adiabatic perturbation in a single-field dominated universe, both in the flat gauge and in 
the uniform energy-density gauge. In order to rewrite the action in a gauge-invariant form 
starting from the uniform energy-density gauge, he needed to introduce a redefinition of the 
first-order perturbations, hence changing their ground state. This redefinition corresponds to 
part of the second-order gauge-invariant curvature perturbation and contributes to the local 
non-Gaussianity. His work was followed by [23] introducing general kinetic terms, [24] for 
two fields and [25-27] for multiple-field models with a generalized kinetic term. In [24] the 
treatment of two fields in the flat gauge showed that no field redefinitions occur (see also [16]). 
Nevertheless, the absence of redeflnitions in this case does not mean the absence of local non- 
Gaussianity, because the action was computed in terms of the scalar fields and not in terms of 
the adiabatic and isocurvature perturbations. This means that in the method of [24] the 6N 
formalism or the long-wavelength formalism is needed to compute the final non-Gaussianity, 
which requires that the slow-roll approximation is imposed at horizon crossing. 

Here we generalize the above results and write the two-field action in terms of the gauge- 
invariant perturbations themselves, both in the uniform energy-density gauge and the flat 
gauge in order to compare the results. We expand the calculation to include second-order 
perturbations and tensor modes, and study the various contributions that occur. Hence 
we derive the exact third-order action, going beyond the slow-roll or the super-horizon ap- 
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proximation. We thus provide the missing tool that will enable people to calculate non- 
Gaussianity, using the in-in formalism [28], beyond these standard approximations used by 
both the long-wavelength formalism and the 6N formalism. 

This paper is organized as follows: in section 2 we provide the gauge-invariant defini- 
tions and conventions for the metric and the field perturbations, along with the description 
of the space-time of the universe, using the ADM formalism. In the first part of section 
3 we study the gauge-invariant curvature and isocurvature perturbations in terms of the 
fields, while next, in subsection 3.2, we make the connection to the fully non-linear spatial 
gradients of the relevant quantities. In the whole of section 3 we use the long-wavelength 
approximation to keep the calculations short and tractable, but we present the generalization 
of the results beyond this approximation in appendix A. In section 4 we construct the exact 
cubic action, going beyond the long-wavelength approximation, to find the redefinitions of 
the perturbations and compare their contributions to the gauge-invariant quantities found in 
subsection 3.1. To keep the main text more accessible, many of the details of the calculations 
have been moved to the appendices. Finally, in section 4.3 we summarize the results from 
section 4, and we conclude in section 5. 



2 Preliminaries 



In this section we give the basic elements required for the calculations in this paper. We 
start by summarizing in subsection 2.1 the ADM formalism, along with the definitions of 
the cosmological quantities, the slow-roll parameters and the field basis we use. In subsec- 
tion 2.2 we provide the conventions of cosmological perturbation theory and clarify different 
approaches in the literature. 

2.1 The ADM formalism 

We will consider a universe filled with two scalar fields with a trivial field metric. The 
generalization to more fields and a non-trivial field metric is conceptually straightforward (see 
[12]), but involves more complicated expressions and calculations. The energy-momentum 
tensor for the two fields (p^ {A, B = 1,2) is 

= <5ABa^</^^9,^^ - g^, Q<5ab<7^"5,¥'^5av^^ + , (2.1) 

where W is the field potential. We will denote the homogeneous part of the fields by (j)^. 
The Einstein summation convention is assumed throughout this paper. We shall work in the 
ADM formalism and write the metric g^i, as 

ds'^ = -N^dt^ + hij{dx' + N'dt){dx^ + N^dt), (2.2) 

where N is the lapse function and A^* the shift. The action takes the form [29] 

S=\j d'^xy/h -2NW+K~^N^\EijE'^-E'^)+NU'^+K~^NR^^^-Nh'^di^'^djipA , (2.3) 

where = SvrG = Svr/mp^, h is the determinant of the space metric hij, R^^^ is the intrinsic 
3-curvature, the tensor Eij (proportional to the extrinsic curvature Kij = —N~^Eij) is 

Eij = I (kj - VtNj - VjN,) (2.4) 
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and n is the length of the canonical momentum of the fields 

= ((^^ - N'di^"^) /N. (2.5) 
Variation of the action with respect to N and A^* gives the energy and momentum constraints 
^-2^(3) -2W- K-'^N-'^{EijE'^ - E"^) - - h'^diif'^djipA = 0, (2.6) 
l{Ei - E6i)] = K^n^^A, (2.7) 



where Vj denotes the covariant derivative with respect to the space metric and E is the trace 
of Eij. 

Following [22] we decompose the space metric as 

/i-^. = a(t)2e2"(*'^)e^'^(*'^\ da'^ = 0, 7^ = 0. (2.8) 

From now on contravariant tensors should be understood as = rf^Ti, where rf^ = 
diag(l, 1, 1), since in the calculations we are showing we have already taken into account 
explicitly the /i*-' part of the initial contravariant tensors. The generalized Hubble parameter 
is defined as 

H=^. (2.9) 

We use the bar for the lapse function, the Hubble parameter and the canonical momentum 
to distinguish these fully non-linear quantities from their background values N(t), H{t) = 
a/{aN) and n(t) = (t)/N, respectively. In this paper we will use as time variable the number 
of e- folds, meaning that d = a, so that the background value of the lapse function N{t) is 
just 1/H{t). 

The background field equation and the background Einstein equations are 

K? K? /H^ \ 
H"^ = -3i?iVn^ - iVW^, H = iVtf, ij2 = _ VW], (2.10) 

with Wa = dW/d(j)^. We construct an orthonormal basis in field space, consisting of 
ef = H'^/n, parallel to the field velocity, and parallel to the part of the field acceleration 
perpendicular to the field velocity [30]. The m = 1 component of physical quantities describes 
the single-field (adiabatic) part, while the m = 2 component captures the multiple-field 
(isocurvature) effects. One can show that for the two-field case the basis vectors are related 
through [15] 

eABcfe^ = -1, (2.11) 

where e^^ is the antisymmetric tensor. The background slow-roll parameters then take the 
form 

e(t) = -^ Jim = £1^51 ^±u) = 2^ ^(t) = }^ + e + J 
dim = li ,^(,\ = ^^^_JL ± 

where Wmn = ^m&n^AB- Throughout this paper the indices m,n will indicate components 
in the basis defined above, taking the values 1 and 2, while i, j are spatial indices and A, B are 
indices of the original fields. In a slow-roll approximation one can think of r]" being related to 
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Wii, T]-^ to W21 and x to W22- The ^ parameters are second-order slow-roll parameters (in a 
slow-roll approximation they are related to the third derivatives of the potential). However, 
we emphasize that we have not made any slow-roll approximations; the above quantities 
should be viewed as short-hand notation and can be large. We also give the time derivatives 
of the background slow-roll parameters and of the unit vectors, 

e = 2iV//e(e + ?7ll), 77" = iVi/ (^^'1 + (t?^)^ + (e - ?7ll)r/ll) , f]^ = N H (^^^ + {e - 27]^^)r]^^ , 
X = NH{ 6r?ll + 2ex - (r?")^ + + + ^r?^^^ + ^ 

= NH { - — + 2r/^e^ + 2eell - 3^" + Oer?" + 3(r?^)" + 3(7?" 





= NH - + 2eC^ - + 9er?^ + 6^^^" " ^^^^ 

et = NHr]^e^, = -NHrj^ef, (2.13) 

where Wi^nn = e^e^efW^ABC- 

2.2 Second-order perturbations and gauge transformations 

In the context of perturbation theory around an homogeneous background any quantity A 
can be decomposed into an homogeneous part and an infinite series of perturbations as 

1 

2' 

where the subscripts in the parentheses denote the order of the perturbation. Up to first 
order the scalar part of the space metric element of (2.8) is equal to 



A{t,x) = Ait) + A^i){t,x) + -A(^2){t,x) + . . . , (2.14) 



hij =a\t){l + 2a^^^)5,j. (2.15) 

When one wants to expand perturbation theory up to second order there are two choices 
found in the literature: either expand (2.15) as Malik and Wands do in [8] to find 

hij = a2(t)(l + 2a(i) + a(2))5ij (2.16) 

or expand directly the space part of (2.8) as for example Lyth and Rodriguez do in [31] to 
find 

h,j = a\t){l + 2a(i) + 2q2^) + a(2))<5ij. (2.17) 

We will take this second approach and use the exponent of the perturbation in our calcula- 
tions. 

Since perturbations depend on the gauge choice we make, we need to construct quantities 
that are invariant under gauge transformations. Under an arbitrary second-order coordinate 
transformation 

the perturbations of a tensor transform as [32] 

^(1) = ^(1) + 

1(2) = i(2) + L^,,,A + Ll^^^A + 2L^(^,i(i), (2.19) 
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where -L/j is the Lie derivative along the vector (3 

{LpAf^l^^l--- = f^^d^A^^ll^l-;: - - • • • + d^.f^^A'^l'^':- + .... (2.20) 

Note here that spatial gradients, having vanishing background values, are automatically 
gauge-invariant at first order, while at second order they transform as 

9il(2) = diA^2) + 2^/3(1) . (2.21) 
3 Super-horizon gauge transformations 

In this section we derive first and second-order super-horizon gauge- invariant combinations. 
We study these during an inflationary period and thus, though we start from the energy- 
density definitions of the perturbations, we naturally end up with field definitions for the 
gauge-invariant perturbations. Our goal is to find the second-order adiabatic and isocurvature 
perturbations in terms of the first-order ones and the slow-roll parameters. 

In this section we restrict ourselves to super-horizon calculations for simplicity, though 
in the next section we will abandon this approximation and study the full action of the 
cosmological perturbations. However in appendix A we present the generalization of the 
results of this section beyond the long-wavelength approximation. We note that the long- 
wavelength (or super-horizon) approximation is equivalent to the zeroth order space gradient 
approximation and is valid once the decaying mode has disappeared (which happens rapidly 
if slow-roll holds during horizon exit), even if there is a subsequent non slow-roll phase. 

In the super-horizon regime one can choose to work in the time-orthogonal gauge where 
N^' = (proof for that choice is given in the next section) and employ the long-wavelength 
approximation to simplify calculations. The latter boils down to ignoring second-order spatial 
derivatives when compared to time derivatives. As a consequence the traceless part of the 
extrinsic curvature quickly decays and can be neglected [33]. Hence the space part of the 
metric can be described by 

/i,,- =a(t)2e2"(*'-)5,,-. (3.1) 

The field and Einstein equations in that case are identical to (2.10), but now the quantities 
involved are fully non-linear. Additionally the momentum constraint (2.7) can be written as 
[12] 

^2 _ 

d,H = -—n^diif^. (3.2) 

3.1 Gauge- invariant quantities 

The well-known first-order adiabatic gauge-invariant curvature perturbation has the form 

NH . , 

Ci(i) = "(I) — ^^(1)' 

where p is the energy density. The subscript without parentheses corresponds to the first 
component in our basis, which is exactly the adiabatic component, while the subscript be- 
tween parentheses denotes the order in the perturbation series. Notice that in the literature it 
is common to work with cosmic time, i.e. N = 1, while the space part of the metric is decom- 
posed using a quantity tp = —a (not to be confused with the ip introduced in appendix C), so 
that the first-order curvature perturbation becomes in that case — Ci(i) = ''/'(i) + {H/ p)p[i)- 
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Here we choose to work with the number of e-folds as time variable so that the first-order 
curvature perturbation is 

Ci(i)=«(i)-^- (3.4) 

The gauge-invariant combination (3.3) is calculated via the requirement that it coincides 
with the curvature perturbation Ci{i) = ck(i) the uniform energy-density gauge where 
P(i) = 0. From now on tilded quantities will denote the uniform energy-density gauge, while 
hatted quantities will denote the flat gauge. One has to use the gauge transformations (2.19) 
for a scalar (here the energy density and the logarithm of the space dependent scale factor 
a) and require that in the uniform energy-density gauge the first-order energy perturbation 
is zero (for details see appendix A). That way one can determine the first-order time shift 
and hence find the gauge-invariant combination corresponding to the curvature perturbation. 
Notice that in the flat gauge, i.e. a(i) = 0, Ci(i) = —P{i)/p- 

Keeping in mind the expansion (2.17), one can repeat the above considerations at second 
order. We find that for super-horizon scales (where we neglect second-order space derivatives 
when compared to second-order time derivatives) the second-order gauge-invariant adiabatic 
perturbation takes the form (see appendix A) 

1. _1. 1 1P(2) , P(i) . 1^(1) P .X 
2^1(2) = r^2) = - 2— + --2 - ■ • (3-5) 

If we chose the second gauge to be flat, i.e. Q(j) = 0, we flnd 

.C.„C,a,.0f. (3.0) 

During inflation we flnd it more useful to work directly with the fields and not their 
energy density, since both the long-wavelength formalism and 6N formalism make use of the 
field values to compute /nl- Using the fields, the first-order adiabatic perturbation becomes 

Ci(i) = 0(1) = a(i) - —eiA^pfi), (3.7) 

since the energy-density constraint /)(i) = is equivalent to eiA'P'^-^ = 0. The detailed 
calculation is shown in appendix B. Notice that in (3.7) we have kept the lapse function N 
arbitrary, as we will also do in all definitions hereafter, but in our calculations H/H is just 
1/0 for the choice N = 1/H. 

The second-order calculation turns out to be more complicated. The details are given 
in appendix B. Here we give the result for the gauge invariant adiabatic perturbation in the 
uniform energy-density gauge and in the flat gauge: 

^Cl{2) = ^"(2) = ^Ql(2) + ('^1(1) " '^2(1)) - ^^Cl(l)C2(l) + Cl(l)Cl(l) 



- d-^d' (C2(1)5.C2(1)) , (3.8) 
where we introduced the auxiliary quantities 

Qm{i) = -r^emAffi) (3.9) 
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and the new combination [12] 



C2(l) = -—e2Affi) = Q2[l), (3.10) 

that represents the isocurvature perturbation to first order. Usually the isocurvature per- 
turbation is described in terms of the gradient of the pressure of the matter content of the 
universe, as for example in [14]. Here we choose to characterize it in terms of the fields them- 
selves and the vector e2A- The latter indicates we are dealing with a purely multiple- field 
effect and hence it is an appropriate quantity to use during the inflationary period to describe 
the non adiabatic perturbations. Starting from the long-wavelength definition of the pressure 
p = — W, one can show that the gradient of the isocurvature perturbation defined in 

[14] is equal to 

r, = d,p - ld,p = -2f/^n2c2i, (s.u) 
p 

where C,2i is the fully non-linear gradient of the isocurvature perturbation (for more details see 
the next subsection 3.2) and f/-*" the fully non-linear generalization of r/-*", i.e. as it is defined 
in (2.12) but with barred quantities [12]. So our definition of C2 agrees with the pressure 
definition of the isocurvature perturbation. The next logical step would be to define the 
second-order isocurvature perturbation C2{2)- However, the above equation shows that there 
is a non-trivial relation between ^2 and the pressure p, involving the non-linear quantities f)^ 
and n, which makes a derivation using the methods of this subsection rather complicated. 
For that reason we prefer to find ^2(2) ™ an easier way in the next section using gradients. 

Notice that unlike in the original definition of Ci(2) in terms of p(2), a non-local term 
appears in (3.8) when one uses the fields instead of the energy density, because of (B.7). The 
time derivatives of the fully non-linear gradients of the perturbations (see next section) were 
found in [12]. Expanding to first order these yield 

Ci(i) = 2r/^C2{i) (3.12) 

for the adiabatic perturbation and 

C2{i) = -XC2{1) (3.13) 

for the isocurvature perturbation, the latter valid only in the slow-roll regime. Then we find 

^~ ^/S I ^ + ^" a2 e-HT?ll -x ^2 , ±/- . /Qi.N 
lf\{2) = 2«(2) = 2^1(2) ^ 2 — ''^(^^ 2 ''i{i)''2(i)' (3.14) 

i.e. without a non-local term. However, we will not use the slow-roll approximation in this 
paper. 

3.2 The gradient of the perturbations 

As an alternative to the Cm defined in the previous section, one can use the gradient quantity 
Q\i along with the isocurvature analogue ^211 both defined in [12] and later in [14] in a covariant 
way, to construct a gauge-invariant quantity. These gradient quantities (not gauge-invariant 
to all orders) are given by 

Cmi = Smidia - r^emAdup'^, (3.15) 
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where now CmA represents the fully non-linear super-horizon version of the orthonormal basis 
vectors, e.g. eiA = n^/n, with IIa = '■Pa/^ since we are working in the super-horizon regime. 
Notice that the basis vectors still obey (2.11) as was shown in [15]. Cmi is by construction 
gauge-invariant at first order, since it has no background value: it is just the gradient of the 
gauge-invariant Cm(i) defined before. 

Expanding to second order we find for the adiabatic perturbation 



1 



1 



1 



^Ch(2) = -^di (a(2) + Qi{2)) - -j^^{i)di^A{i) - 2-7eiA9?(i)9j<3i(i)- 



(3.16) 



In the uniform energy-density gauge this gives (see appendix B) 

2^11(2) = 2'^«"(2)' 

while in the fiat gauge where diOt = and Ci{i) = Qi{i)^ we find 



(3.17) 



Cli(2) 



di 



Qi 



+ 



e + 7]» 



■Ci(i)C^jCi(i)i 



Ci(i) - C2(i) j - ??"^Ci(i)C2(i) + Ci(i)Ci(i) 



C2{l)9iC2{l) 

(3.18) 



where we used the basis completeness relation and (B.8) to rewrite the terms. Cii(2) iii 
the uniform energy-density gauge (3.17) coincides with the gradient of the gauge-invariant 
second-order adiabatic perturbation. However, by comparing (3.8) to (3.18) we see that 
in the fiat gauge Cii(2) is the gradient of the gauge-invariant curvature perturbation (1^2) 
expressed in the fiat gauge plus a new non-local term. This is in agreement with the findings 
in [14]. This new term is nothing else but the gauge transformation of Cii{2)- A quantity 
with zero background value as Q is, transforms as (2.21). One can check, using the gauge 
transformations (2.19) for p and for a and requiring that /0(i) = 0, that starting from a fiat 
gauge and transforming to the uniform energy-density gauge the time shift is T(i) = (^1(1) 
(see appendix A), so that 

2Cmi{2) = 2^™(2) + Cl(l)Cmj(l)- (3.19) 

Next we try to find the second-order gauge-invariant part of the isocurvature perturba- 
tion by expanding (3.15), 



^C2i(2) = ■7^diQ2(2) + -J^^^^'^B{l)di^A{l) " '^^(ilA'^fl)diC,2(l), 



(3.20) 



or equivalently. 



^C2i(2) = l^diQ2{2) - -je2A^fi)diQi{i) - ^eiA^fi)diC2{l), 



(3.21) 



where we used (2.11) to express 62 in terms of ei. The first-order uniform energy constraint 
diQi{i) = alone implies 



7)^2H2) — di 



2^2(2) - ^^2(1) 



+ Cl(l)t?jC2(l), 



(3.22) 
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where we used (B.14) for the two last terms in (3.21). For the flat gauge diCt = we find 
using (B.8) 




Cl(l)'9iC2(l)- 



(3.23) 



We notice that the term in the second fine corresponds again to a gauge transformation 
famihar from the curvature perturbation case studied earher (3.19). In appendix A we verify 
that the rest of the expression is a gauge-invariant quantity corresponding to the one in 
(3.22). Indeed this expression is gauge-invariant beyond the long- wavelength approximation 
as shown in appendix A. 

We conclude that the gradients of the perturbations are in some sense equivalent to the 
perturbations themselves, since both allow for the definition of gauge-invariant second-order 
adiabatic and isocurvature quantities. However, since the gradients are defined using fully 
non- linear quantities, their equations of motion can be treated more easily, as was shown in 



4 The cubic action 

An alternative way to calculate the second-order gauge-invariant quantities and reconsider 
their meaning, is to compute the third-order action for the interacting fields. Maldacena 
[22] was the first to perform that calculation for a single field, in the uniform energy-density 
gauge. In this way he managed to find the cubic interaction terms due to non-linearities of 
the Einstein action as well as of the field potential, which among other consequences change 
the ground state of the adiabatic perturbation Ci(i)- This change can be quantified through 
a redefinition of the form [22] 



where Cic is the redefined perturbation. One sees that the correction term of the redefinition 
coincides with the surviving quadratic term of the single-field limit of the transformation 
(3.8), taking into account that the super-horizon adiabatic perturbation is constant in that 
case. In [22] the curvature perturbation was considered a first-order quantity, while the 
second-order curvature perturbation was not taken into account, since its contribution in the 
uniform energy-density gauge is trivial: it introduces a redefinition of the form Ci(i)+Ci(2)/2 = 
Cic(i) (foi' proof, see subsection 4.2). Seery and Lidsey [24] performed the same calculation 
for the multiple-field case in the fiat gauge in terms of the scalar fields (p^ and not of the 
adiabatic and isocurvature perturbations Cm- They found no redefinitions, but as mentioned 
before their results would have to be supplemented by the 6N formalism (with its associated 
slow-roll approximation at horizon-crossing) to say anything about the non-Gaussianity of 
the gauge-invariant perturbations Cm- 

In this section we generalize the above calculations to second order in the expansion 
of the curvature perturbation in both the uniform energy-density gauge and the fiat gauge. 
Doing so we compute the full form of the third-order action. The latter not only consists 
of the cubic interactions of the first-order curvature perturbations, but also of lower order 
interaction terms of the second order quantities. We first perform the calculation relevant 



[12]. 



Ci(i) — Cic(i) + 



2 
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to the first-order quantities and then add the second-order effects. In this section we only 
present the scalar part of the action, but in appendices C and D the tensor part can be 
found as well. We emphasize that in this section we no longer make the long-wavelength 
approximation, so that the results are valid at any scale. 

4.1 The second-order action 

We start by performing our calculation in the gauge eiA'P'^^ = 0. This constraint reduces 
to the uniform energy-density gauge outside the horizon, which is why we will continue to 
refer to the tilded gauge as the uniform energy-density gauge. Prom now on we drop the 
explicit subscript (1) on first-order quantities. We will keep this part brief since its results 
are already known, but we give the basic elements of the calculation in appendix C. The 
second-order action takes the form 



^2 = Jd^xL2 

= jd^x e|- a-| ({dCif + {dC2f) + a'H(cl + C| - 4r/^CiC2 + 2XC2C2) (4.2) 

+"'^(/?^ - 2e^ - {vh' + Hv^f + - 3e(ryll - x) + 2r/llx)cl|, 

where L2 is the second-order Lagrangian. While we have started from an action describing 
the evolution of the fields d and e2A'f'^ we have now constructed an action in terms of the 
adiabatic and isocurvature perturbations and (2- The equations of motion that Ci and (^2 
obey are {5L2/6C,rn being a short-hand notation for the relevant variations of the Lagrangian) 



Ci + ( 3 + e + 277" ) Ci - 2r/^C2 + 2(-C^ - 2e77^ - 377^)C2 + 2a-dXi 



-^{2a^Hd'^X) + 2a-^52Ci = 0, (4.3) 
d.t H 



(>C2 



C2 + ( 3 + e + 2??ll ) C2 + 2r/^Ci + (e" + 2e^ + AerjW + 3x)C2 +2a-9^C2 = 



e 



where d^X = e^i — 2erj-^C,2 (for the reason of introducing A see appendix C). Thus we have 
found the evolution equations for the first-order adiabatic and isocurvature perturbations. 
Their super-horizon limit coincides with the equations derived in [12] for the gradient of the 
perturbations, since up to first order Cmj(i) = f^iCm(i)- One can show that the first-order 
energy constraint, which outside the horizon reduces to 

Ci - 2??^C2 = 0, (4.4) 

is the first integral of the super- horizon part of the first equation of (4.3), i.e. without the 
space gradient. In fact it was shown in [34] that this is the case at all orders. In the same 
paper it was found that assuming the slow-roll limit, ^2 = — XC2 is the super-horizon first 
integral of the equation for (^2, which can be easily verified. 

While working in the flat gauge we find the same action (4.2) (see appendix C). So 
the curvature perturbations Cm satisfy to first order the same equations in both gauges as 
expected, due to the gauge invariance of Cm (or equivalently the gauge invariance of the 
action). 
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4.2 The third-order action 

In this section we compute the third-order action. Again we present only the final results, 
while in appendix D we give the intermediate steps of the calculation. In the same appendix 
we also give the tensor-scalar part of the action. The scalar cubic action in the uniform 
energy-density gauge due to the first-order perturbations Cm takes the form 

-53(1) = 'S'3{i) - / d^a;^/m (4.5) 

with 

/i = '-^C! - v^CiC2 + CiCi - 4^ ((5Ci)' - d-^d^d^d,Cid,Ci)) 

/2 = (e + r/")CiC2 + C2Ci + ^C?. (4.6) 

The exact form of 5'3(i) can be found in appendix D or equivalently it is the cubic part of 
(4.15). The reason for introducing 5'3(i) without the tilde will become clear below. 

The terms proportional to 5L2/6(^m, i.e- the first-order equations of motion, can be 
removed by a redefinition of Cm [22] and lead to a change in the ground state of the pertur- 
bations. This works as follows. The cubic terms of the action (i.e. 5'3(i)) are not affected 
by the redefinition, because the redefinition always involves terms proportional to C,^, which 
would give quartic and not cubic corrections. It is only the second-order terms (i.e. 52) 
that change. Indeed one can show that under a redefinition of the form = Cmc + fm, the 
second-order action changes as ^2 = 52c + {5L2/5C,m)fm- These new terms cancel out the 
relevant terms coming from the cubic action (remember that the total action up to cubic 
order is the sum of the second and third-order action) and we are left with 

5*3(1) = 5'3(i)(Cmc). (4.7) 

If we repeat the same calculations for the flat gauge (see appendix D), performing several 
integrations by part, we find that 

^3(1) =53(l)(U). (4.8) 

This is a consequence of the action staying invariant under a gauge transformation. Never- 
theless if one associates the redefinition appearing in the uniform energy-density gauge to a 
change in the ground state of Cmj it would mean that directly after horizon crossing, when 
super-horizon effects have not yet been switched on, the second-order contribution to C,i would 
be zero for the flat gauge and non-zero for the uniform energy-density gauge. In terms of 
non-Gaussianity, this can be restated as: the non-Gaussianity present after horizon-crossing 
is different for the two gauges. Indeed if one was to calculate the three-point functions for 
the above action, one would need to perform two steps. First, change to the interaction 
picture, where it can be proved that the interaction Hamiltonian up to and including cubic 
order is just Hint = —Lint-, where Lint are the cubic terms of the Lagrangian, and compute 
the expectation value (CcCcCc) as in [22]. Second, take into account that the flelds have been 
redeflned as = Cc + • Then the three-point correlation function can be written as 



(CCC) = (CcCcCc) + 2A[(CcCc)(CcCc) + cyclic]. 



(4.9) 




These new terms, products of the second-order correlation functions, are only present in the 
uniform energy-density gauge if we restrict ourselves to £'3(1) . 

In order to cure this bad behavior we need to add to the above results the effect of the 
second-order fields. We find (see appendix D) 

= I d^-^ ^ ( + ^ ) + ^%^}- (4-10) 

Since all terms in 53(2) are proportional to 5L2/6Cm, S^2) only contains redefinitions of Cm- 
Notice that the second-order lapse and shift functions do not appear in the final action, since 
these two are multiplied by a factor equal to the energy and momentum constraint equations 
(C.2). On the other hand, the second-order field perturbations are dynamical variables that 
obey second-order equations of motion that cannot be set to zero in the action. The single- 
field limit of this action is just the term proportional to Ci{2)) since Qm{i) = identically 
in that case for the uniform energy-density gauge. The term proportional to Ci(2) i^i 5*3(2), 
along with the terms proportional to A in 5*3(1) , originate from the contribution of A^* in the 
action. The latter vanish outside the horizon since then d^X coincides with the super- horizon 
energy constraint and hence is identically zero. So if we were to study only the quadratic 
contributions of the first-order perturbations outside the horizon, we would be allowed not 
only to ignore the tensor parts of the metric [33] , but also work in the time-orthogonal gauge 
N' = 0. 

Coming back to the redefinition, its final form, including the tensor parts (see appendix 
D), is 

C2 = C2c - ^ + C2C1 + ^C? + {e + r/")CiC2. (4.11) 
Finally we perform the above calculations for the fiat gauge and find the action 

The redefinitions in the flat gauge take the simple form 

Ql{2) 



Ci = C: 



Ic 



2 



C2 = C2c-^. (4.13) 

We want to write these redefinitions as well as the action itself in terms of gauge-invariant 
quantities and compare them. We would also like to compare with the definitions of the 
second-order gauge-invariant perturbations found in appendix A and section 3. After using 
the second-order uniform energy constraint (B.13) and the uniform energy gauge definition 
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of (^2(2) (3.22) we can rewrite (4.11) as 



Ci + ^ = Cic + CiCi + ^ (C? - Cl) - v\iC2 - d~'d^ (C23.C2) - \d~\%d^dKi) 

C2 + ^ = C2c + C2C1 + C2C1 + - '^2) + + ^")ClC2 - 5-25* . (4.14) 

When comparing the first equation of (4.14) with (A. 14), we see that we recover (4.13). The 
same is true for the isocurvature part of the redefinition: comparing the second equation of 
(4.14) with (3.23), we recover the redefinition for C,2 (4.13). Hence the two redefinitions are 
the same, as is necessary for the action to be gauge-invariant. Notice that the single-field 
limit of (4.14) is Ci + Ci(2)/2 = Cic + (e + r/ll)(^^/2 in agreement with the total redefinition 
found in the uniform energy-density gauge. 

Equation (4.14) is the implicit definition of the redefined, gauge-invariant C,rnc- One 
can see that up to and including second order, it is a function of only the combination 
Cm.{i) + Cm(2)/2. One can also notice that the purely second-order perturbation Cm(2) does 
not occur explicitly in the cubic action (see e.g. (4.15) below). Hence once could in principle 
consider the quantities Cm{2) (and similarly Qm(2)) as auxiliary quantities and try to avoid 
introducing them in the first place, but consider the quadratic first-order terms directly as a 
correction to the first-order perturbations, as is done for the single-field case in [16]. While 
the calculations would be roughly equivalent, we have chosen not to follow this route for two 
reasons. In the first place it seems conceptually simpler to us to expand the perturbations 
and the action consistently up to the required order, and more logical to view quadratic first- 
order terms as a correction to a second-order quantity than to a first-order one. Secondly, 
in the multiple-field case (as opposed to the single-field case), one would have to introduce 
the second-order quantities at some intermediate steps anyway in order to find the correct 
non-linear relation between the Qm and Cm (which is derived from the second-order gauge 
transformation performed in appendix A). 

So in the end we have managed to find the source of the non-Gaussianity present at 
horizon crossing due to first-order perturbations and identify it with the quadratic terms 
of (4.14). With source here we mean the second-order perturbation that, when contracted 
with two first-order perturbations, gives the bispectrum. The super-horizon limit of (4.14) 
was derived and used in our previous paper [15], but here we have not only generalized the 
result, but also have obtained a much better understanding of the gauge issues. Equation 
(4.14) is gauge- invariant, as it should be. Additionally, the redefinition of the perturbations 
that we perform is essential not only to simplify calculations but also to find the gauge- 
invariant form of the action itself. We clearly see that the quadratic corrections in the flat 
gauge seem to be zero if one takes into account only the first-order fields. In that gauge 
all of the second-order contributions are hidden in the second-order fields as opposed to the 
uniform energy-density gauge where part of the quadratic contributions is attributed to the 
redefinition of the first-order Cm and the rest of them lie in the second-order field. 

4.3 Summary 

Let us summarize the results of this section. Cosmological gauge-invariant perturbations 
should obey a gauge-invariant action. Using first-order perturbations the action up to third 
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order is the same in the uniform energy-density gauge and the flat gauge only after a re- 
definition of Cm in the uniform energy-density gauge Cm = Cma + /ml (4.6) (the subscript 1 
indicating the use of only first-order perturbations) and takes the form 



S — S{Crn) — 'S'(Cmci) 

(1 + eCl) (C! + Cl) - ( C2diC2 + Ci^.Ci 



+2C2 (xC2 + eCi ((e + ,")C2 + V^Ci)) + Cl (-4??^C2 + eC?(3r/" + 2e) 




X 



/2eW22i 




+ e -3(?7^)2 + (e + ?7ll)2 +3e(x-e-r/ll) Ci 



(4.15) 



where we have kept the notation d^X = e^i — 2er/^C2 in order to mark clearly the terms that 
vanish outside the horizon, namely the terms proportional to A along with the terms involving 
second-order space derivatives. This is one of our main results. We managed to compute 
the cubic action for adiabatic and isocurvature perturbations in the exact theory, beyond 
any super-horizon or slow-roll approximation. Its single-field limit coincides with the action 
computed by Maldacena in [22] or by Rigopoulos in [16]. Let us examine the implications 
of this action. Forgetting about the redefinition of the perturbations in the uniform energy- 
density gauge, the form of the action is gauge- invariant. One can use it to easily calculate the 
non-Gaussianity related to the interaction terms as is explained in detail in [28, 35]. This is 

(3) 

known in the literature as fj^^, the parameter of non-Gaussianity related to the three-point 
correlation function of three first-order perturbations, which is only non-zero in the case of 
intrinsic non-Gaussianity. 

However, taking into account the need for a redefinition of the perturbations in the uni- 
form energy-density gauge, one might worry that the action is not actually gauge-invariant. 
The action in the uniform energy-density gauge before the redefinition has extra terms that 
are proportional to the second-order equations that the perturbations obey. This means that 
when calculating the non-Gaussianity in the uniform energy-density gauge, one not only has 
contributions due to the the interaction terms in the cubic action, but also ones due to the 
redefinition of Cm? which contribute as explained in (4.9). They are part of what is known in 
the literature as /^|^, the parameter of non-Gaussianity related to the three-point correlation 
function of a second-order perturbation (in terms of products of first-order ones) and two 
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first-order perturbations, which reduces to products of two-point functions of the first-order 
perturbations. 

This would mean that the non-Gaussianity calculated in the two gauges would not be the 
same due to the lack of any redefinition in the flat gauge. However, if one takes into account 
only the corrections coming from first-order perturbations, the redefinition associated to the 
second-order perturbation is not complete as one can check by comparing the super-horizon 
version of the adiabatic part of (4.6) with (3.8). As we showed, the solution of this issue is 
to include second-order fields since they also contribute to the cubic action. As one would 
expect these do not change the action itself, so that (4.15) still holds. The effect of the new 
terms is to redefine the perturbations in both gauges. It should be noted that, if one had 
incorporated all quadratic first-order terms (found by a second-order gauge transformation 
as in appendix A) directly as a correction to the first-order perturbations, one would have 
found the two contributions 5*3(1) and 53(2) together and hence there would have been no 
initial discrepancy between the two gauges. However, we explained at the end of section 4.2 
our reasons for proceeding in this way. So in any case we finally obtain 



where (mc is given in (4.14). Now the two redefinitions as well as the action itself are the 
same for the two gauges, hence the action is truly gauge-invariant and the fj^\, related to 
the products of first-order Cm in the redefinitions, is the same in the two gauges. 

This exact action allows one to compute without the need for the slow-roll approx- 
imation at horizon crossing that is essential for both the long-wavelength formalism and the 
5N formalism: the long-wavelength formalism needs slow-roll at horizon crossing in order to 
allow for the decaying mode to vanish rapidly, while the 6N formalism requires it in order 
to ignore the derivatives with respect to the canonical momentum. Additionally, up to now 
only the slow-roll field action [24] (and not the action of the Cm themselves) was known, so in 
order to compute the non-Gaussianity at horizon crossing one had to use the long-wavelength 
or 5N formalism to transform to Cm and hence one was in any case required to make the 
assumption of slow-roll, even if the exact action for the fields would have been known. It will 
be interesting to investigate models that do not satisfy the conditions for the long-wavelength 
or 5N formalism using the action (4.15). 

In order to connect the redefinitions to some previously derived results in the literature 
we assume the super-horizon and slow-roll approximations. The super-horizon approximation 
is already assumed in (4.14) and it can be supplemented by the condition Ci = 2r]-^(^2- The 
slow-roll assumption translates into ^2 = — xC2- Then the quadratic part of the redefinitions, 
relevant to fi^l, takes the form 



The redefinitions in this form were used in [15] to find the second-order source term of the 
evolution equations for the super-horizon perturbations. Their contribution to the super- 
horizon /^l^ was calculated in that paper using the long-wavelength formalism. In the equal- 
momenta case it was shown to be 




(4.16) 




(4.17) 




1 + {vi2y 



(4.18) 
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where the index * indicates the time when the scale exits the horizon and vi2 is essentiahy a 
transfer function showing how the isocurvature mode (2 sources the adiabatic mode Ci (see 
[15] for details, where this term is part of what is called gsr)- Directly after horizon crossing 
or equivalently in the single- field limit, when V12 = 0, this reduces to the well-known result 
by Maldacena — = e* + 



5 Conclusions 

In this paper we settled some unresolved issues concerning gauge invariance at second order 
in inflation with more than one field. Although the gauge-invariant curvature perturbation 
defined through the energy density has been known for many years, the energy density is 
not the quantity that is used in calculations of inflationary non-Gaussianity. These use the 
scalar fields present during infiation instead of their energy. We found this gauge-invariant 
quantity in terms of the fields and discovered that it contains a non-local term unless slow-roll 
is assumed. 

We have also managed to make contact between gauge transformations and the redefi- 
nitions of the curvature and isocurvature perturbations occurring in the third-order action. 
Since [22] it has been known that the redefinition of the curvature perturbation in the action, 
introduced to remove terms proportional to the first-order equations of motion, corresponds 
to its gauge transformation. However, these terms appear at first sight to be absent in the 
flat gauge which would have had as a consequence the absence of quadratic contributions 
of first-order curvature perturbations at horizon crossing in this gauge and hence a gauge 
dependence of the related horizon-crossing non-Gaussianity (using Wick's theorem one can 
calculate the three-point correlation function due to these terms, as we did in [15]). We 
have extended the calculation for both gauges to second order and proved that in both of 
them the contributions are the same. The difference is that, in our perturbative approach, 
in the uniform energy-density gauge a part of these contributions is due to the first-order 
corrections and the other part to the second-order fields while in the flat gauge they are all 
due to the second-order fields. 

In addition to the adiabatic one, we also found the gauge-invariant isocurvature pertur- 
bation defined in terms of the scalar fields by studying the relevant fully non-linear spatial 
gradient defined in [12]. Usually isocurvature perturbations are studied in terms of the 
pressure of the fields. Following [12] we found a definition using the fields themselves that 
demonstrates the orthogonality of this quantity to the curvature perturbation. While rewrit- 
ing the action, these isocurvature perturbations appear naturally in the form we have defined 
them, thus showing that this quantity is the relevant one to use during inflation. 

In order to achieve the above we computed in section 4 the exact cubic action for the 
perturbations, going beyond the slow-roll or super-horizon approximations (in appendix D we 
also give the tensor part of the action). This can prove very useful for future calculations. Up 
to now one had to impose the slow-roll condition at horizon crossing in order to calculate the 
non-Gaussianity. This was because the only two-field action available was that of the fields 
given in [24], thus demanding slow-roll at horizon crossing in order to be able to use the long- 
wavelength formalism or the 6N formalism to find the curvature perturbation bispectrum. 
The action we provide here can be used directly with the in-in formalism [28] in order to 
calculate the exact non-Gaussianity beyond any restrictions, slow-roll or super-horizon. 
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A Gauge transformations 



From the infinite number of possible gauge-invariant combinations, we choose to work with 
quantities constructed from the energy density and the logarithm of the space dependent 
scale factor a. We will consider a gauge transformation /3(j) = (r(j),^(j)) from the hatted 
gauge to the tilded gauge, where for the moment both gauges are taken to be arbitrary 
(not yet the flat and uniform energy density gauge). Notice though, that the space part of 
the transformation is not relevant outside the horizon, since when introduced in the rela- 
tions below, it is connected to a second-order space derivative [8]. Within the super-horizon 
approximation, we find using (2.19) 

P{1) = P(l) + P^(\)^ P{2) = P{2) + PT{2) + T(l) (2^(1) + pT(i) + pT(i)^ , 

"(1) = "(1) "(2) = "(2) +^(2) +^(1) (2a{i) . (A.l) 



We want to construct a gauge-invariant quantity that reduces to a(j) in the uniform energy- 
density gauge, which we now identify with the tilded gauge so that yO(j) = 0. This way we 
find 

-'(1) = -'(2) = — 1 {^■^) 

and obtain 

Ci(i) = a(i) = "(1) - — , 

V - ^/^ 1P(2) , P{1)^ 1^(1) P ..ON 

Z ^ ' Z ^ ' Z ^ ' Z p P ^ P P 

Notice that the initial hatted gauge is still arbitrary, but if one was to associate it with the 
fiat gauge a(j) = 0, then the time shift would become T(i) = Ci(i)- 

Next, we derive the exact gauge-invariant adiabatic and isocurvature perturbations, 
going beyond the super-horizon approximation. We use (2.19) for the scalar fields and the 
space part of the metric tensor (2.8), and find 

'^fi)=^fi)+'^%) (A.4) 
+ \\m = + + ^(i)'^ii + ^i^i^Ci) + \ + • (A.5) 

Here we followed [32] and split the component i of the space shift as qi = dig + q;^, where 
d^q^- = 0. We choose the uniform energy-density gauge, defined by eiyit^^^ = and use (A.4) 

to find the first-order time shift to be T(i) = — eiA<^(^)/<^. Then the trace and the traceless 
part of (A.5) give 



"(1) = «{!) + ^(1) + (A-6) 



+ 2 (9,9, - i<5,,-52) + d,qfi^^ + d.qf,^^. (A.7) 



In order to make the definition of the super-horizon adiabatic perturbation at first order 
(A.l) to agree with (A. 6), we choose = (any choice of ^(i) is a gauge-invariant quantity. 
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but only q^^^ = corresponds with the adiabatic perturbation Ci{i) i^i t^is hterature). Note 
that while working with the super-horizon approximation, no such choice needs to be made, 
and remains arbitrary in that case. Similarly, we also assume that g^^. = 0, so we find 

Ci(i) = = a(i) - ieiA<^fi) and 7(i)ij = 7(i)ii = 7(i)jj- (A.8) 

Using (A. 4), one easily finds that the isocurvature perturbation at first order is gauge- 
invariant since e2A4>^ = 

C2(i) = -ie2A'^fi) = - Je2A'^fi)- (A.9) 

We now fix the hatted gauge to be the fiat one, a(j) = 0, in order to lighten the calculations. 
This implies that T(i) = Ci(i)- 
At second order we find 

^(2) = "^m + ^'^^ + + Ci(i)0^ + 2^fi)) , (A. 10) 

where we have omitted a term proportional to which as mentioned above is chosen to 

be zero. At second order we choose the gauge ^Qi(2) = ^^^2~^2{i) ~^ d~'^d'^ ^C2(i)t^iC2(i)^ > see 
(B.13), that reduces to the uniform energy-density gauge on super-horizon scales (for the 
definition of Qi{2) in (B.ll) and details about that gauge choice, the reader can refer to 
appendix B). Using (B.8) we find from (A.IO) 

^(2) = Qi(2) + (e+^") (Ci(i) - C2(i))+Ci{i)Ci(i)-2??^Ci{i)C2(i)-29-'a^ (C2(i)9^C2(i)) • (A.ll) 

Before turning to the adiabatic perturbation, let us prove that the first line of (3.23) 
is a gauge-invariant quantity corresponding to the one in (3.22). This is true in the exact 
theory, beyond the long-wavelength approximation, for g^x) = 0. Multiplying (A.IO) with 
— e2A/(2</'), noticing that e2A4>^ = Oi and using (B.8), one finds 

1 ~ 1 - T]^ 

2<32(2) = 2'32(2) + ^Ci(i) + C2(i)Ci(i) + (e + ??")Cl(l)C2(l) (A.12) 

and indeed by comparing the total gradient of (3.22) and (3.23) we see that it corresponds 
to the second-order gauge-invariant isocurvature perturbation. 

For the second-order adiabatic perturbation we need to perform the gauge transforma- 
tion (2.19) of the space part of the metric tensor between the uniform energy-density gauge 
and the fiat gauge 

"{2)'^ij + \l{2)ij = \l{2)ij + T(2)Sij + didjq{2) + ^ (di(lt2)j + 9jqf2)?) + Cl(l)7(l)ij 

+diC,i(i)dj\ + 9jCi(i)<9iA + Ci(i)Ci(i)5ii - ^^2^^»'^i(i)^i'^i(i)' (^-13) 

where we substituted T^j^) = C,i(^i-^ and T(^2) is given in (A.ll). In order to find Ci(2) we subtract 
the trace and the d~'^d^d^ of (A. 13) to eliminate the term proportional to ^(2) and obtain 
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the result 

2Ci(2) - 2«(2) 

= 2'3l{2) + Cl(l)Cl{l) + (''1(1) ~ '^2(1)) - ??^Cl(l)C2(l) - d~'^d'' (C2(l)5iC2(l)) 

+ \{d'\diQi - d-^d'd^idiXdjCi))- (A.14) 

This is the second-order adiabatic gauge-invariant perturbation in the exact theory, the 
generalization of (3.8). 



B Super-horizon calculations 

In the first part of this appendix we give the detailed calculations of subsection 3.1, while in 
the second part we show the ones relevant to subsection 3.2. 

An important property of the long-wavelength assumption is that outside the horizon 
the uniform energy density can be recast in terms of the fields at least at first order: one can 
show that the exact Oi-Einstein equation (3.2) outside the horizon can be rewritten as 

dip = -2,HIiBdi^^ . (B.l) 

Again p denotes the fully non-linear energy density p = n^/2 + W . Expanding (B.l) to first 
order and using the background equations to prove that 

p = -m^NH, (B.2) 

one can show that outside the horizon 

M = — ei^<^^ (B.3) 

p UN ^(1) ^ ' 

so that 

Ci(i) = 0(1) = a(i) - — eiA(/?fi). (B.4) 

and thus the energy-density constraint p(i) = is equivalent to eiA'^^^-^ = 0. 

Unfortunately the nice property described by (B.3) does not hold anymore at second 
order. After expanding up to second order, combining the second-order Einstein equations 
and using the completeness relation of the field basis, one can show that (note that the zeroth 
order lapse function is taken from now on to be N{t) = 1/H(t)) 

^-f - -^e,A^t2) = {j'^Mi)) - ^) + (B.5) 

with 

]^diA = 'n-^e2B<ffi)eiA0i(pfi) + e2B0i(pfif2A'pfi), (B.6) 

where we used (2.13) to simplify the expressions in terms of the slow-roll parameters. We 
see that the purely second-order contribution of p(2) is recast as a second-order contribution 
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of <^^2)^ some quadratic first-order terms and a non-local term arising essentially from the 
Oi-Einstein equation. In the flat gauge diOi = this non-local term can be written as 



'n Ci{i)C2{i) + ^ C2(i) 



+ C2(l)(5iC2(l)i 



(B.7) 



where when needed we employed the following useful relations (valid for the flat gauge beyond 
the long- wavelength approximation): 



1 



-eiA<^(i) = Ci{i) + (e + ??")Ci(i) - C2(2), 

-ie2A'^fi) = C2{i) + (e + ??")C2(1) + ??"^Cl(2), 
derived by differentiating the adiabatic perturbation and the new combination 



C; 



H 



2(1) 



(B.9) 



that represents the isocurvature perturbation to first order. Putting everything together in 
(3.6) we find the second-order gauge-invariant curvature perturbation in the flat gauge to be 



Ql{2) + Cl(l) 



TjCl(2) = TJ"(2) 



- d-'d' (C2(1)5.C2(1) j , 

where we defined the auxiliary quantities 

Qm{i) = 



e + ^2 ±A /- I / / 

— ^ — C2(i) - V Ci(i)C2(i) + Ci(i)Ci(i) 



(B.IO) 



(B.ll) 



We turn now to the calculations relevant to the gradient of the perturbations. In the 
uniform energy-density gauge we can use (B.l) to find the constraints 



1 



3 







1) 

P{2) = -^4'Adi(ff2) - ^'fA{l)di'P(l) 



or diQ^i) = 



or 



^diQn2) - -J^'^A{l)di(pfi) - —^Y~diC2(l) + C2(l)'9iC2(l)> 



(B.12) 



(B.13) 



so that when inspecting (3.16) we see that 

To derive the last equality in (B.13) we used the completeness relation of the field basis along 
with the following relations valid beyond the long-wavelength approximation: 

1 



;eiA<^fl) = -??"^C2(1), 



1 



-762^(^(1) = C2(i) + (e + ^")C2(i), 



(B.14) 



derived by differentiating the definition of the isocurvature perturbation and the first-order 
uniform energy-density gauge constraint eiA'P'^-^ = 0. 
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C Second-order action calculation 



In order to rewrite the action we first need to calculate the extrinsic curvature. To do that 
we decompose N = 1/H + Ni, = d^ip + N]_^ where diN\_ = 0. From now on we drop 
the explicit subscript (1) on first-order quantities and set = 1 to lighten the notation 
(notice though that the final results remain unchanged when we restore = SvrG, since all 
are absorbed in e when rewriting the fields in terms of C,m)- We start by performing the 
calculation in the gauge 

eiA^^ = 0, (C.l) 

which we call uniform energy-density gauge, since the above constraint reduces to zero energy 
perturbation outside the horizon. We first use the energy and momentum constraint (2.6), 
(2.7) to find that to first order 



^ H if' 



0, 



(C.2) 



It turns out that we do not need to calculate the shift or the lapse function to higher order, 
since in the action those terms are multiplied by constraint relations and hence vanish. 

We start by working out the scalar part of the action. Keeping in mind the gauge 
constraint (C.l) we perturb (2.3) to second order 



S2 = ^ I d'^xla^^" 



A,~B 



+H{1 - HNi + H'^Nt - H 



N2 

2 



—ae 



+ Ni)2 {{da f + 2d^a) + ^di^Ad'^^ 



H 



6(1 + d)2 + + 2<j)^^A + 

(C.3) 



where we have omitted a total derivative with respect to tp. We then use the background 
Einstein and field equations to eliminate some terms and find that the term proportional to 
A''2 vanishes. Now the second-order action can be written as 



^2 = ^/d^x<; a-'H 



H 



Ci(^-4^<^A + </''Cij +9Ci'(-6 + 
2e{dCif + d'^Adi^^'^ 



1 



Wab(p^^^ - 36CiCi + 



(C.4) 



The terms of (C.4) proportional to can be recast in terms of the curvature perturbations 
by applying the completeness property of the field basis and (B.14), so that after integrating 
by parts and using H = —eH it can be written as 



So 



d^x e< a^H 



Ci + Cl - VC1C2 - 3(x - e - ryll)C| + 2(e + ?7ll)C2C2 



2/-2 



1 r 



(aci)2 + (8(2)' 



(C.5) 
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or after further integration by parts 



S2 



d^x e -al ( {ddr + {dC2f ) + a'H[c:i + (i - 4r?^CiC2 + 2XC2C2 



|^-2.^-(r?l')^+3(r?^)2 + ^r?^e^-36(r?ll-x)+2r?llx)c||(C.6) 



We can reach the same result while working in the flat gauge diOt = 0. One can prove 
that Ni = —eQi/H, N]_ = and d'^'ip = d'^X = eC,i — 2er]^(2- The ^jJ terms cancel out and the 
second-order action takes the form 



-a^d'i^Adi^' 



(C.7) 



Using the definition of Cmi along with the background equations (2.10), (2.12) and (B.8) this 
can be rewritten as (C.6). 

The second-order tensor part of the action in both gauges takes the form 



5- 



27 



d'^i27 = \ I d'x{^H{^,,f - ^(9,7.,)'}, 



(C.8) 



where L27 is the second-order Lagrangian for the tensor modes. We also give the equation 
of motion of the gravitational waves 

SL2-y 



5% 



(0,9) 



which we are going to use in the next section. In this paper we will not discuss the evolution 
and physics of gravitational waves, but at linear order this is a standard subject in the 
literature, for a discussion see for example [29]. 



D Third-order action calculation 

In order to compute ^3 we follow the same procedure starting from the uniform energy- 
density gauge. Notice that -/V3 will multiply {—2W + 6H^ — iP) in exact analogy with iV2 
in S2, so it vanishes. Moreover, the overall factor multiplying A'^2 is the first-order energy 
constraint (C.2), so it can be consistently set to zero as well. 

We start by computing the cubic action of the first-order curvature perturbations up to 
A^i involving only scalar quantities 



1 



'53(i) = 7T ld''x{a-'e 



(— + Ni 



2W - 2Wa^'' - Wab^^^^ - ^Wabc^^^^^^ 



1, 



+H 



(1 - HNi + H'^N't - H'-^N^) ( - 6(1 + q)^ + + 2<j)'^^A + ^ 



,~2 



+ 



-ae 



(d'd^didji: - {d^i>f^ (1 - HNi) - Ad'ipdiCid'^ip - 2i^Ad'i>diip^ 

(D.l 



[j^ + iVi) ( a'^A^^^^ + + 2{dCif 
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After using the background equations and the definitions of the perturbations, eq. (D.l) 
takes the form 



5. 



3(1) 



Jd''xf^'e'^H[e{l - Cl) (C? + Cl + 2(e + rj\\)U2 + {{v^f + {e + r?")') Cl - 27/^C2Ci) 



-2e7?^C2Ci - 3(1 + CMx - e - 7?II)C| + 



+i (d'd^ijd.djij - {d^f) (1 - Cl) - 2eav((e + 77ll)C25iC2 + C25iC2 



1 



-a-(Ci + Cl) 29^C + {dCiY + e('9C2)^ 



(D.2) 



By performing integrations by parts in (D.2) we find 



<i(C? + Cl) - 2Ci5U9,Ci - 2C25Ua,C2 - 2[er,\\ + {r,\\f + {^^f)CiU2 
+ (3ery^ + e)ClC2 + 2(67?^ + ^^)CiC2Ci - (e" + 2er,\\ + {r,\\f + {v^f^ClCi 



+ (2e^ + 3e,ll -{r^Wf -{.^f + e")c?Ci + fi'^Cl 
+ ^2(e + 7?" )e^ + r?^ (2e(3 + e) + 677" - - 3x) ^ C?C2 

+ ^- e(4e2 + 6e + 12er/ll + ??" (9 + Sr/") + 8(??^)2 + 2^" - ^x)+^ 



2e W221 



-3(77" r - ^iv^r - 2?7"C" - 2?7^C^ CiCi - 2(e + v'K2d'm2 



,2\\2\ 



+ ae^^Ci 



{dCif+{dC2f 



-f (^C? - .-C1C2 + CiCi - ^(aCi)^ + ^^-^^^^(^.Cia.C.) 
+^5Xi9.A - ^d-^d'cP{d,XdjCi)) - ^ ((^ + '?")CiC2 + C2C1 + ^Ci')(k3) 



where 6L2/SCm ai'e the first-order equations of motion. We can further integrate by parts the 
rest of the action to simplify it and prove that it takes the form of the flat gauge action (D.7), 
as expected since the action should be gauge- invariant. The terms involving A along with the 
terms with space gradients vanish outside the horizon in the long-wavelength approximation, 
since A is equal to the first-order super- horizon energy constraint (4.4). 
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Finally we include the second-order fields. The extra terms in the action are 



^3(2) = 2 Id^xWe''^^ 



{j,+N,){-WA^^^)-WABV>f2)^'') + H{l-HN^){cPA^f^)+^A^f2)) 



H 



(D.4) 



where ip^ without a subscript always denotes the first-order perturbation. After performing 
integrations by parts we find 



c - I A^^j '^-^2 ( Ci(2) Qi(2)\ SL2 Qm_ \ 

^3(2)- + ^ 



K2 2 ]• 



(D.5) 



Next, we perform the same calculation for the flat gauge, starting from 



53(1) = - U^'xla^ 



+Ni)(^-2W- 2Wa(p^ - Wab^^^"" - ^WABC^^if''^^) (D.6) 



+H{1 - HNi + H^Nf - H^Nf) (^-6 + <j)^ + 2<j)^i^A + ^"^ + 49^^ 



again taking into account that N2 multiplies the first-order energy constraint and thus we 
set it to zero. We find using the definition of Cmj along with (2.10), (2.12) and (B.8) 



S- 



3(1) 



<i{Ci + C2) - K2d'Xd,C2 - 2Qid'\diCi 
+2e{e + r?ll)CiC2C2 + 2er/^CiC2 + 6(37?" + 2e)C?Ci 



+ 



27^221 



+ e (-3{ri^f + (e + rj^^f) + 3e{x - e - r?ll) ) (ICi 



-2(e + Tj^%d'Xd^C2 + 4r?^C23Ua,Ci + ^Ci{d' XdiOjX - {d^Xf 
+^Ci((5Ci)^ + (5C2)^)|. 



(D.7) 
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Finally we include the second-order fields. The surviving terms in the action are 



53(2) = § Id^xla^ 



H 



-H^'^d'i;di(pA{2) 



a—d'ip di(pA{2) 



(Di 



and they can be rewritten as 



'3(2) 



4 j 6L2 Ql{2) SL2 Q2(2) 



(D.9) 



In the last part of this appendix we consider the tensor scalar part of the action. There 
will be no contributions from the second-order fields, since these cancel due to 'jij being 
transverse. We start from the action for two scalar and one tensor modes in the uniform 
energy-density gauge 



d^x 



H L 



' - (3Ci - Qi)iijd'&ii^ + dk^ijd'd^^d^] }, 

which after integrations by parts becomes 

1 



(D.IO) 



'CC7^ 



+ 



H 

6L2J 



a 



dXid^Ci + dX2d'C2 (D.ll) 



In the flat gauge one can find directly after substitution in (2.3) the first line of (D.ll), so 
that there are no redefinitions. 

Finally we calculate the part of the action consisting of one scalar and two tensor modes, 
starting from the uniform energy-density gauge 



5, 



1 



C77 



a 

AH 



iCi + Ci)idkli,f}iB.12) 



or equivalently 



'C77 



6L 



27 



;Cli^^,f - \%d,y^d'x] + ^<,(4^,,.)')P-13) 

while in the flat gauge we find directly 

^C77 = IJ d^x{a3/? [iCi(%f - l%dkY^d'x\ + -^eCiid^j,,)']. (D.14) 

The three tensor modes action does not contain any redefinitions. For details the reader 
may look in [22]. 
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